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Abstract. Using STIS on board the HST we have obtained a spectroscopic map of the bipolar jet from RW Aur 
with the slit parallel to the jet axis and moved across the jet in steps of (/.'OT. After applying a velocity correction 
due to uneven slit illumination we find signatures of rotation within the first 300 AU of the jet (I'.'S at the distance 
of RW Aur). Both lobes rotate in the same direction (i.e. with different helicities), with toroidal velocities in the 
range 5-30 km s~^ at 20 and 30 AU from the symmetry axis in the blueshifted and redshifted lobes, respectively. 
The sense of rotation is anti-clockwise looking from the tip of the blue lobe (P.A. 130° north to east) down to 
the star. Rotation is more evident in the [OI] and [Nil] lines and at the largest sampled distance from the axis. 
These results are consistent with other STIS observations carried out with the slit perpendicular to the jet axis, 
and with theoretical simulations. 

Using current magneto-hydrodynamic models for the launch of the jets, we find that the mass ejected in the 
observed part of the outflow is accelerated from a region in the disk within about 0.5 AU from the star for the 
blue lobe, and within 1.6 AU from the star for the red lobe. Using also previous results we estimate upper and 
lower limits for the angular momentum transport rate of the jet. We find that this can be a large fraction (two 
thirds or more) of the estimated rate transported through the relevant portion of the disk. The magnetic lever 
arm (defined as the ratio rA/ro between the Alfven and footpoint radii) is in the range 3.5 - 4.6 (with an accuracy 
of 20 - 25 %), or, alternatively, the ejection index ^ = dln(Macc)/dr is in the range 0.025 - 0.046 (with similar 
uncertainties). The derived values are in the range predicted by the models, but they also suggest that some 
heating must be provided at the base of the fiow. 

Finally, using the general disk wind theory we derive the ratio / Bp of the toroidal and poloidal components of 
the magnetic field at the observed location (i.e. about 80 - 100 AU above the disk). We find this quantity to be 
3.8 ± 1.1 at 30 AU from the axis in the red lobe and -8.9 ± 2.7 at 20 AU from the axis in the blue lobe (assuming 
cylindrical coordinates centred on the star and with positive z along the blue lobe). The toroidal component 
appears to be dominant, which would be consistent with magnetic coUimation of the jet. The field appears to be 
more tightly wrapped on the blue side. 

Key words. ISM: Herbig-Haro objects — ISM: jets and outflows — stars: formation — stars: pre-main sequence 
— stars: individual: RW Aur 



1. Introduction 

The collimated Herbig-Haro (HH) jets observed on parsec- 
scale lengths in star formation regions are always associ- 
ated with young stellar objects (YSOs) that are still in 
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their accretion phase. Accretion and ejection of matter 
are believed to be intimately related phenomena, through 
the presence of a magnetized accretion disk. The complex 
interplay between accretion and ejection is modelled using 
several theoretical approaches (cf. Camenzind et al. ll990[ 
Ferreira lTMTl Konigl & Pudritz 2000. Shu et al HUM and 
references therein) . These models have in common the idea 
that the jet is generated through the interaction of plasma 
with rotating magnetic field lines that are anchored to the 
star/disk system. The fluid particles lifted from the disk 
are forced to slide along the rotating magnetic field lines, 
and are accelerated and collimated into bipolar jets. The 
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most attractive aspect of this approach is the fact that 
the magneto-centrifugal scenario at the same time justi- 
fies the acceleration of the jets and the extraction of the 
excess angular momentum from the disk. This mechanism, 
in combination with disk viscosity (the nature of which, 
however, is still uncertain), should contribute to slow the 
disk material down to sub-Keplerian rotation. In this way 
matter is allowed to move radially toward the central star 
and finally accrete onto it. 

This compelling theoretical picture, however, has re- 
ceived little in the way of observational confirmation, as 
pointed out by Eisloffel et al. f2000|). The main reason is 
that the process occurs at small distances, corresponding, 
for the nearest star formation regions (130 - 150 pc), to 
subarcsecond scales. 

The "core" of the central engine lies within a few AU 
from the source, which is not spatially resolvable with cur- 
rent observational techniques. Important constraints on 
the launching region, however, can be deduced from the 
quantities observed in the first 100-200 AU of the flow, in 
which the jet achieves its collimation and its final poloidal 
velocity. This region has been accessed recently with high 
angular resolution observations of jets from evolved T 
Tauri stars. Such objects have been observed from space 
with the Hubble Space Telescope Imaging Spectrograph 
(HST/STIS, on resolution, e.g. Bacciotti et al. 12013 
Woitas et al. l2U(J2|l . and from the ground using large tele- 
scopes with Adaptive Optics (e.g. Dougados et al. .2000(1 . 

Among the properties investigated in the above stud- 
ies, a special place is taken by the detection and analy- 
sis of the rotation of the jets around their symmetry axis 
(Bacciotti et al. Testi et al. IMHl Anderson et al. 

l^nnSl Coffey et aL OUMl Pesenti et al."2DU3). These stud- 
ies have shown that the measured toroidal velocities are 
similar to the values predicted by the disk-wind magneto- 
centrifugal models. The obtained measures of rotation, 
however, also help to constrain the properties of the ac- 
cretion/ejection "machine". For example, in combination 
with mass density estimates from line ratios, rotation pro- 
files can yield the amount of angular momentum carried 
away from the accreting system by the jet. This can then 
be compared with the disk accretion properties. As an- 
other example, the application of special conservation laws 
can in principle give information on the properties of the 
magnetic field in the jet, a quantity that is fundamental to 
all jet acceleration models, but that is notoriously difficult 
to directly examine by observation. 

The aim of this paper is to illustrate the potential of- 
fered by such a combined observational/theoretical inves- 
tigation to constrain the accretion/ejection structure at 
the origin of the jet. First, we present results from a new 
observational study of the RW Aur jet, conducted with 
multiple slits oriented along the flow axis (Sections 2 and 
3). The adopted technique is similar to the one used in our 
flrst rotation study, that concerned the jet from DG Tauri 
(Bacciotti et al. l2002|) . Here, however, the technique is ap- 
plied for the flrst time to a bipolar jet, on a more extended 
and more finely sampled region. In addition, more details 
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Fig. 1. Schematic drawing of the observing mode for the 
bi-polar jet from RW Aur. The STIS slit is kept parallel to 
the flow axis, and stepped sequentially in seven different 
positions (labelled 51, S2....S7) each time by O'.'OT. The 
slit width is O'.'l, so the apertures are slightly overlapping. 
The sense of rotation derived in the present paper is illus- 
trated by the arrows around the symmetry axis. 

are given about the correction routines for uneven slit illu- 
mination. Then, in the second part of the paper (Section 
4), we combine the obtained results with morphological 
and excitation properties of this jet derived on subarc- 
second scales from previous studies and partly based on 
the same dataset (Woitas et al. 120021 , hereafter Paper 

1, Dougados et al. 12000(1 . In this way we show how, us- 
ing general theoretical principles, one can derive physical 
quantities that are crucial for a description of the accre- 
tion/ejection region close to the origin of the jet. Finally, 
in Section 5 we summarize our conclusions. 

2. Observations and data analysis 

The HST/STIS observations of RW Aur were carried out 
on 10 December 2000. Details of the data acquisition and 
reduction were described in Paper I. Briefly, the dataset 
consists of seven spectra (hereafter 5*1,..., 57) taken at slit 
positions parallel to the jet axis (P. A. — 130°, Dougados 
et al. l2000|l . This observing procedure is shown in Fig.^ 
The pixel scale is O'.'05/pixel for the spatial axis and 
0.554 A/pixel for the dispersion axis. The latter value cor- 



Woitas et al.: Jet rotation and derived physical parameters in the RW Aur jet 



3 



responds to « 25 kms ^ /pixel. The slit width is O'.'l. The 
covered spectral range contains Ha and the prominent for- 
bidden emission Hnes [OI]AA6300,6363, [NII]AA6548,6583 
and [SII]AA6716,6731. However, blueshifted [OI]A6300 
emission cannot be used for the present analysis because 
the line profile is cut by the detector edge. 
In the course of the observations the slit is moved trans- 
versely across the jet in steps of 0'.'07 from southwest {SI) 
to northeast (S'7). The position of 5*4 corresponds to the 
jet axis. If the jet rotates, one expects non-zero peak ve- 
locity differences S'7 — 5*1, 5*6 — S2 and 5*5 — S3 that have 
the same sign at all separations from the star. However, 
one has to take into account that uneven slit illumination 
may introduce spurious velocity shifts causing peak ve- 
locity differences similar to those described above, which 
may then be present also in the case of a non-rotating jet. 
This effect has been eliminated applying correction rou- 
tines developed by one of us (A.M.), which also have been 
used to study the rotational motions in the DG Tau jet. 
The calculations performed in the correction routine are 
illustrated in Appendix A of Bacciotti et al. (|2002) and 
in Marconi et al. H2003P , and the procedure can be briefly 
summarized in three steps as follows. 

(i) In the first step a bidimensional surface Ix{x,y) is 
calculated that simulates the observed brightness distri- 
bution in each emission line at the positions (x,y) on the 
sky intercepted by the seven STIS slits. This is done by 
fitting separately for each line a series of 2-D surfaces with 
an elliptical base simulating the knots along the flow to 
the line fluxes measured along the seven slits at once. The 
best fit is found through an iterative algorithm. It is then 
necessary to assume a model radial velocity field of the 
nebula u{x,y). In our case we wish to determine the spu- 
rious velocity caused by STIS, so we assume an arbitrary 
constant value u, equal for all the slits and positions on the 
sky, i.e. a non-rotating uniformly moving jet is assumed. 

(ii) In the second step the routines calculate the av- 
erage velocity measured at the detector v\{x,y) for each 
slit, line and position considered. This is done by convolv- 
ing the model velocity field as it would be seen at the 
detector without optical distortions with the model of the 
brightness distribution obtained at step (i), with the Point 
Spread Function of HST at the wavelength of interest and 
with the slit aperture offset with respect to the axis of 
the jet as in the observations. Details of this calculation 
can be found in the Appendix of Bacciotti et al. (|2()t)2(l . 
Since the model velocity of the gas is set to be the same 
value for all the slits, the searched instrumental spurious 
velocity is Vx{x,y) — v\{x,y) — u, the result being inde- 
pendent of the chosen value of the model radial velocity. 
The absolute values otV\{x,y) range between 2 and 8 km 
s~^ and have opposite sign depending on the relative po- 
sition of the slit and the jet axis, since slits to the left or 
to the right of the axis have opposite illumination gradi- 
ents. The values for Vx{x, y) are higher (in absolute value) 
where the illumination gradient is steeper across the jet, 
i.e. typically at the position of the knot peaks. 



(iii) Finally, the spurious "shift" is calculated between 
the considered opposed slits as, e.g. V\{S1, y) — V\{S1, y). 
The calculated shifts are instrumental offsets, and there- 
fore they are subtracted from the raw shift measurements, 
to cancel the effect of the uneven slit illumination. 

The spurious velocity depends only on the gradient of 
the slit illumination, and not on the integrated flux across 
the slit, so any systematic error in the photometry would 
not affect our determination. Therefore, the uncertainty 
associated with the correction depends only on the accu- 
racy with which the shape of the observed emission dis- 
tribution is fitted by the assumed 2-D surfaces convolved 
with the instrumental response. The other parameters en- 
tering the calculation (see Bacciotti et al. I2()()2|l . i.e. the 
instrument PSF, the detector scale factor and the line flux 
measured by STIS are all known with a much higher ac- 
curacy than that associated with the fitting procedure. 

The accuracy of the correction has been estimated 
a posteriori by investigating the variations in the deter- 
mined spurious shifts due to an imposed change in the 
scale factor of the brightness distribution of the surface 
fitting each knot. Variations of the adopted illumination 
gradient producing a change by 10 - 15 % in the value of 
the fitted fiux, which is about the accuracy of our itera- 
tive fitting algorithm, lead to a change of 5 - 8% in the 
determination of the spurious velocity shift. Thus we take 
this factor as the uncertainty of the correction, which, in 
the case of the largest spurious shifts 16 km s~^), cor- 
responds to less than ±1.5 km s^^. Also, we checked that 
imposed small misplacements of the slit positions or incli- 
nation (up to 10%) do not lead to any appreciable change 
in the determination of the spurious velocity. 

It should be noted that in the real case in which the ve- 
locity field is not constant, an additional deformation is in- 
troduced by the passage of the light through the telescope 
and spectrograph. In our case, however, this contribution 
can be neglected because we are actually interested in spu- 
rious velocity differences between symmetrically opposed 
slit pairs, and the additional contribution would cancel in 
the difference because the absolute value of the real ra- 
dial velocity, as well as the intensity, have a symmetrical 
variation with respect to the jet axis or central slit. 

In practice, the spurious velocity is positive for the 
slits SI, S2 and S3, located south-west of the axis, and 
negative for 5*5, S6, S7. So for the adopted configuration 
(S'7 — SI, S6 — S2, S5 — S3) the instrument produces a 
spurious negative shift in both jet lobes. The effect is more 
pronounced for the forbidden lines (FELs) than it is for 
Ha. This is due to the fact that the bright Ha emission is 
more evenly distributed across the transverse direction in 
the area covered by the seven slits (see Paper I). Similarly, 
other detected differences in the spurious shifts for the 
different lines can be attributed to the different transverse 
spatial distribution of the emission (see Table 2 and its 
description in Section 3). 

After subtracting the spurious negative shifts the ro- 
tation signature becomes more evident for the outermost 
slit pairs (larger positive shift), while in the innermost 
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Fig. 2. Example of a Gaussian fit to the line profile that 
yields the peak velocity. This plot shows [SIIJA6731 at slit 
position S3 and at a separation of (/.'Q from the star. 

pair the measured raw negative shifts are turned back to 
weakly positive or null (i.e. consistent with zero inside the 
measurement error of ± 7 km s^^). In some positions of 
the innermost pairs, however, the net shift in [SII] and Ha 
lines remains significantly negative. For an explanation see 
Section 3. 

The validity of the illumination correction can be 
tested using the Ha emission at the stellar position. As 
this emission is close to saturation, the line profiles will 
be dominated by the HST PSF and are not supposed to 
show any rotation signatures. Any velocity shifts between 
lateral slits will be due to instrumental effects. The find- 
ing that there are no significant velocity shifts in Ha close 
to the star after applying the illumination correction (see 
Fig-EJ thus strongly indicates that this correction routine 
really removes velocity shifts caused by uneven slit illumi- 
nation. 

The peak velocities of individual lines in all spectra and 
in both outflow lobes were determined from Gaussian fits 
to the line profiles. Many T Tauri jets show two veloc- 
ity components of forbidden line emission, a high velocity 
component (HVC) with Vr ~ 100 — 200kms~^ and an 
additional low velocity component (LVC) with a typical 
radial velocity of 5 to 20kms"i (Hirth et al HW/jl . The 
HVC is thought to come from a coUimated jet close to the 
star, while the LVC might be the signature of a poorly 
coUimated disk wind (Kwan & Tademaru 1 1 9881 [TMK| . In 
Paper I we have, however, demonstrated that there is no 
significant separate LVC in our data. Therefore fitting sin- 
gle Gaussians is sufficient. As an example we show in Fig.|5] 
the fit to the [SHJA6731 line at slit position S3 and sepa- 
ration 0'.'9 in the redshifted lobe. The uncertainty of the 
peak velocities is about ±5kms^^, and the uncertainties 
of velocity differences are thus ±7kms~^. 

For some selected positions and emission lines we have 
derived velocity dilferences also by means of a cross- 
correlation technique that directly measures the displace- 
ment of two line profiles from different slits, with a typ- 



ical accuracy of 5kms~^ on the velocity difference. This 
method is in principle more robust than Gaussian fitting 
as no special shape of the line profile has to be assumed. 
Its usefulness is limited, however, in situations where the 
velocity range is only marginally sampled, and this is the 
case here with the narrow emission lines of the RW Aur jet. 
Thus cross-correlation does not lead to a distinct improve- 
ment, but the results are consistent with those obtained 
from the Gaussian fits. Reliable measurements are not pos- 
sible for all slit positions and separations from the source. 
The RW Aur jet has a FWHM of about 20 AU within 
its first arcsecond and is thus more strongly coUimated 
than the jet of DC Tau (FWHM « 50 AU at d = O'.'S, see 
Paper I, Fig. 4 therein). For this reason, there are cases 
where there is no significant signal at distinct separations 
in 5*1 and 5*7. In [Nil], which traces emission very close to 
the jet axis, this will even be the case for more central slit 
positions. Furthermore, quenching effects make the FELs 
extremely weak within O'.'l - 0'.'2 from the origin, thus we 
exclude this region from the analysis. Although the jet is 
clearly visible up to « 4" (« 800 AU) from the star in 
the STIS data (Paper I) we restrict the rotation analysis 
to the first « 300 AU approximately. Over this separation 
range our data will reflect signatures from the initial jet 
channel, whereas further away interaction between the jet 
and the circumstellar environment might become impor- 
tant. To obtain de-projected separations we assume the 
inclination angle of the RW Aur jet to be i = 46 ±3° with 
respect to the line of sight (Lopez-Martin et al. I2003|l . In 
Paper I we adopted i — 53°, but this estimate was based 
on the proper motion of only one jet knot and furthermore 
affected by an erroneous pixel scale given by Dougados et 
al. f20001 for their 1997 observations (the corrected value 
is given by Lopez-Martin et al. i2003) . With the new value 
of the inclination angle, 1'.'5 on the sky corresponds to 
300 AU along the outflow direction. 

3. Indications of rotation 

Fig.O shows the peak velocities for all slit positions and 
at different separations from the star. These velocities are 
corrected for uneven slit illumination (see Sect.|2J) and 
for the heliocentric radial velocity of RW Aur A, i. e. 
-t-23.5kms^^ as derived from the Li I fine at A « 6707A 
from our data. We note that Petrov et al. H2001|l repeat- 
edly observed RW Aur using high resolution echelle spec- 
troscopy and found a mean heliocentric u,. of +16kms^^ 
with periodic variations of ± 6 kms~^ (P « 2.8 d). Such 
variation in the adopted systemic velocity does not af- 
fect our determinations of velocity differences between the 
two sides of the flow. For the present analysis, however, 
in which we also give transverse proflles of the absolute 
radial velocity, we decided to adopt the former value for 
self-consistency. 

The resulting velocity differences S'7 — 51, 56 — 52 and 
55 — 53 are plotted in Figs .0] and as a function of the 
separation from the star. The result is not as impressive 
as for DG Tau, due to the fact that the RW Aur jet is 
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Fig. 3. Peak velocities for different slit positions and separations from the star. Crosses denote [SIIJA6731, squares 
[SII]A6716, circles [OI]A6363, triangles [OI]A6300, plus signs [NII]A6548, asterisks [NII]A6583, and diamonds Ha. 
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Fig. 4. Peak radial velocity differences in the redshifted 
outflow lobe as a function of separation from the source. 
Upper panel: difference between the values measured at 
slits 57 and S'l, or at 30 AU from the axis. Middle panel: 
5*6 - S2, corresponding to 20 AU from the axis. Lower 
panel: S5 - S3, or 10 AU from the axis. The plot symbols 
have the same meaning as in Fig.|3| The large diamonds 
represent Ha at the stellar position, the open one before 
the correction for uneven slit illumination and the filled 
one afterwards. 



more collimated (see below). One can see, however, that 
in all slit pairs there is a definite tendency for the velocity 
shifts to be located on the same side of the zero line in 
both lobes. The scatter of the data points in the different 
lines is sometimes large, and it is probably due to different 
excitation conditions, as discussed later. Nevertheless, the 
averages of all velocity shifts in all lines and separations 
from the star show the same sign in all slit pairs and in 
both lobes (Table 1). Thus, the emission from the north- 
eastern parts of both lobes is more redshifted than the 
emission from the south-western parts. This leads us to 
suggest that both outflow lobes rotate clockwise looking 
from the tip of the redshifted lobe down to the star. The 
matter thus flows in the jet with opposite helicities. These 
results are fully consistent with the study in Coffey et al. 
(2004). 



20 



ot +/- 0. '21 
from jet oxis 



°* ^ o 
o o 

o □ ^ o 

o o o 



-10 



20 



ot +/- 0."14 
from jet oxis 



X o o n X * ^ 

■ ■■0 



xo 00 



X xo «H o 



20 



r 



-10 



at +/- 0. '07 
from jet oxis 



o^oO 



* + ° M o 

V X X <>o 

X ^ «53o X o 00X0 o 

it X M 0*0 Ox*<>S^-n 

O w 



0.0 



0.2 



0.4 0.6 0.8 1.0 1.2 
Arcseconds from source (blue lobe) 



1.4 



Fig. 5. Same as Fig.0I but for the blueshifted outflow 
lobe. 

Table 1. Mean velocity differences (in kms~^), averaged 
over all lines and all separations from the star 





57 -SI 


56 




52 


55 


- 53 


Redshifted lobe 


12.9 ± 1.7 


0.2 


± 


0.8 


2.0 


±0.6 


Blueshifted lobe 


5.2 ± 1.3 


1.1 


± 


1.6 


3.8 


±0.9 



Given the fact that the measured rotation velocities 
are of the order of only w lOkms^"'^ one has to look care- 
fully at artifacts that might mimic rotation. The problem 
of uneven slit illumination has already been discussed in 
Sect.|21 Another issue that has to be taken into account 
here is that the slit orientation might not be exactly par- 
allel to the jet axis. Gaussians fits to the emission profiles 
transverse to the jet, as well as the output of the illumina- 
tion correction routines, do indeed show that the slit set is 
well-centred on the star, but the chosen slit P.A. (130 deg) 
turns out to be larger by 1-2 degrees than the real P.A. of 
the jet direction (blue lobe - the same misalignment is also 
noted in Coffey et al. 120041 ). This causes the slits of one 
given "pair" to trace gas with different poloidal velocity, 
which will lead to apparent velocity shifts even if the jet is 
non-rotating. From Fig. 3 one can estimate that the cor- 
responding "false" shifts would be 1 - 5 kms^^, which is 
not negligible. For the given position angle the mimicked 
rotation would, however, have the opposite direction in 
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Table 2. Velocity shifts (in kms in different emission 
lines, averaged over all separations from the star. 
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7.5 


1.6 


[SII]A6716 


4.1 


8.1 


2.6 


[SII]A6731 


-5.6 


-6.4 


3.3 



comparison to the observed motion. This means that the 
true rotation velocities are in fact somewhat higher than 
those suggested by Figs. 0] and |5l and given in Tabled 
strengthening our result even further. 
Finally, the observed velocity shifts could in principle re- 
flect a situation where the poloidal velocity field is not 
symmetric to the jet axis, but without any rotation. This 
can happen in asymmetric bow shock wings, but is very 
unlikely to occur along a distance of 300 AU, which is 
much larger than the knots of the RW Aur small-scale jet 
(Paper I). 

As mentioned above, the data show a notable scatter in 
Figs. 4 and 5, and there are several possible explanations 
for this, dumpiness, spatially unresolved shocks and vari- 
ations of the poloidal velocity may play a role in produc- 
ing the observed scatter. Some scatter would be present, 
however, even if the flow was laminar with the poloidal 
velocity smoothly decreasing from the axis to the jet bor- 
ders ( "onion- like" kinematic structure), as predicted by 
the MHD stationary models. According to these models 
the layers located closer to the axis rotate faster, so dif- 
ferent emission lines that trace different regions along the 
line of sight will also show different radial velocity shifts 
if measured at the same position on the sky. From Table 
12 one can see that the rotation signatures are more pro- 
nounced in the [01] and [Nil] lines than in [SII] and Ha. 
We have indeed found that the transverse FWHM of the 
emission in [01] and [Nil] is smaller than the FWHM in 
[SII] (Fig. 4 of Paper I) . The emission in Ha also presents 
a larger transverse extension. Thus at the same position 
across the jet, corresponding to two coupled slits, the "[01] 
jet" and the "[Nil] jet" are observed closer to their borders 
than the [SII] (and Ha) jet. Pesenti et al. I|2()()4(l note that 
projection of the emission along a line of sight that crosses 
the beam of the jet tends to cancel out the radial velocity 
differences produced by rotation, this effect being more 
severe for distances closer to the axis and lower spatial 
resolution of the instrument. Therefore, at any point of 



observation, the rotation signatures will be more evident 
in [01] and [Nil] , as the emission is less averaged across the 
line of sight in this case. The same projection effects can 
explain the fact that the evidence for rotation is weaker 
in the innermost slit pairs than it is in 57 — 51. Note 
also that some values are negative in Table 2 for [SII] and 
Ha. We interpret this as follows: We have seen above that 
the set of slits is subject to an unwanted inclination of 
about 2 degrees with respect to the jet axis, in such a way 
that the observation of a non-rotating jet would produce a 
'false' negative shift. On the other hand, we have pointed 
out above that for any given slit pair the "[01] jet" will 
show larger rotation shifts than the "[SII] jet", because 
of projection effects combined with a different spatial dis- 
tribution of the emission. As explained above, in [01] the 
rotation signature dominates over this inclination effect. 
In [SII] lines, instead, the rotation signature is very weak 
at the same position. This is because the aforementioned 
enhanced projection effects for the [SII] lines overwhelm 
the rotation here, and the net effect is the observation 
of a negative shift in [SII]. Ha may also be affected in 
the innermost slits. We have not attempted to correct for 
this effect, however, since the value of the "false" rota- 
tion depends on the kinematic model assumed for the jet. 
Pesenti et al. (|2()()4jl also predict rotational velocity shifts 
for "cold" and "warm" disk wind solutions. Our observed 
rotation velocities are intermediate between these two pre- 
dictions, but closer to the warm disk wind solution. This 
point will be discussed further in Section 4. 

Finally, we note that a decrease of the velocity shift 
with distance from the star is marginally seen in some of 
the lines and slit separations, as, for example, in the red- 
shifted lobe for [OI]A6300 in the pair 57 - 51 and 56 - 
52, and for [SII]A6716 in the 56 - 52 pair, while for the 
blueshifted lobe a decrease is apparent in the Ha points for 
the 57 - 51 pair and for [SII]A6716 in 56 - 52. We note 
that such a behaviour is predicted by all the models of 
magneto-centrifugal acceleration as a consequence of the 
conservation of kinetic angular momentum at large dis- 
tance from the Alfven surface (see below) and the widen- 
ing of the flow. Given the scatter of our datapoints and 
their expected accuracy, however, we refrain from present- 
ing a quantitative analysis of this aspect. Instead, in the 
rest of the paper we will use the radial velocity shifts mea- 
sured in the first 0'.'6 - O'.'S from the star to derive as much 
information as possible about the physical properties of 
the launching region of the jet. 

4. Properties of the accretion/ejection structure 

From our results one can derive quantities useful to con- 
strain the properties of the jet acceleration region. In par- 
ticular it can be checked if the values obtained are consis- 
tent with the self-similar disk-wind models described, e. g. 
in Konigl & Pudritz 1,2000.) and Ferreira ^UU^ . 

The launching region of the jet is sketched in Fig. 
This figure illustrates the configuration of the nested mag- 
netic surfaces attached to the inner portion of the accre- 
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tion disk. The star is at the origin of a cylindrical co- 
ordinate system (r, z,0), and the flow is assumed to be 
steady-state, axisymmetric and to satisfy the ideal MHD 
equations. The poloidal component of the magnetic field 
Bp = BrBr + Bze^ is described by Bp = Va xe^/r, where 
a{r,z) = Const labels the magnetic surfaces, i.e. the sur- 
faces enclosing a constant magnetic flux. The poloidal ve- 
locity and magnetic field are easily shown to be parallel, 
and they are related by the expression: pVp = A; (a) Bp, 
with k constant along each magnetic surface. So the mag- 
netic and flow surfaces are coincident and can also be la- 
belled with their "footpoint" radius r = ro, or the distance 
from the star on the disk where the surface is anchored. 
The mass in the visible jet is ejected from a region of the 
disk ranging from an inner radius ro,in to some outer ra- 
dius rg oxt to be determined. The relevant mass fluxes are 
denoted by large grey arrows in the figure. 

The disk and the rigidly anchored field lines rotate 
rapidly, and the centrifugal force makes the fluid parcels 
lifted from the disk surface (by the thermal pressure) be 
flung out along the open field lines. The matter flowing 
in the jet trails behind the fleld on a given surface. This 
generates a toroidal field component and, as a con- 
sequence, a "magnetic torque" that brakes the disk, and 
extracts energy and angular momentum from it. In the 
acceleration process, the matter reaches a point where 
the magnitude of the poloidal velocity equals the Alfven 
poloidal velocity, or \vp\ = ca = \Bp\/y/iTTp [p is the mass 
density). The loci of such points constitute the so-called 
Alfven surface. In the self-similar disk wind models (see, 
e.g., Casse & Ferreira l2000|l . this surface is conical, and its 
section is also indicated in Fig.Elwith a dashed line. Above 
the Alfven surface, the inertia of the matter overcomes the 
magnetic forces, and the intensity of the toroidal fleld is 
increased substantially. This in turn generates a magnetic 
force ("hoop stress") directed toward the axis, that can 
coUimate the flow. 

The Alfven surface lies at a few AU above the disk, 
a region that is still not resolved by observations. The 
datapoints we obtain in the first 0'.'2 - O'.'S from the star, 
however, correspond to the region immediately above, and 
can be used to investigate the properties of the gas that 
has just been accelerated and coUimated. The steps of 
the adopted procedure are illustrated in what follows. 
Measured and derived quantities are summarized in Tables 
El and respectively. 

4.1. Footpoint radii 

The inner edge of the ejection region in the disk is com- 
monly identified with the point where the disk is trun- 
cated by the interaction with the stellar magnetosphere. 
The corresponding disk annulus rotates at the same rate 
as the central object, and we will take this "corotation ra- 
dius" Tc as a fiducial value for ro,in- For a typical T Tauri 
star rc ~ 0.03 AU (Shu et al. .2000)1 . 




Fig. 6. Sketch (not to scale) of the configuration of the 
flow/magnetic surfaces (solid thin lines) attached to the 
inner section of the accretion disk around RW Aur A, 
according to the disk- wind scenario. Thick grey arrows 
indicate the mass flows. Note the different distances of the 
footpoints for the blue and red flows seen at 0'.'4 above the 
disk. This is due to the asymmetry in the bipolar jet. 



An estimate of the outer radius ro,oxt of the ejection 
region for the matter seen in optical lines can be derived 
from the combination of the poloidal and toroidal veloci- 
ties {v-p, v^) we measured in the external layer of the flow 
located at robs from the jet axis. In a magneto-centrifugal 
wind the footpoint radius of the jet component located at 
Tobs from the axis with measured and Vp can be esti- 
mated with the relationship valid at large distances from 
the source provided by Anderson et al. H2003|l : 

— (lStj)"(loi^)" 

VlOOkms-i/ \1 MqJ ' ^ ' 

in which we will take M^, = 1.3 ± 0.1 Mq for the mass 
of RW Aur A (Woitas et al. I2001|l . This simple equation 
is valid under the condition <C Wp, which in the con- 
text of disk winds is achieved as soon as r > va, where 
the 'Alfven radius' va is the distance from the axis at 
which the magnetic/flow surface rooted at tq intersects 
the Alfven surface. We note that at the sampled distances 
above the disk (50 - 100 AU), this asymptotic regime may 
not have yet been entirely reached by the outer stream- 
lines of the flow. It can be shown, however, that in the 
case of 'light' TTS winds, like the one treated here, the 
results given by the simplified formula above do not dif- 
fer by more than a few percent from those obtained by 
a rigorous treatment of the conservation equations (see 
Anderson et al. H2003|l ). 

According to the expression above, in the red lobe, and 
for Tobs = 30 AU (slits S7—S1), one has Wp^Rcd 130 ± 10 



Woitas et al.: Jet rotation and derived physical parameters in the RW Aur jet 



km s^^ (Paper I). For the toroidal velocity in the first 0'.'2 
- O'.'S from the source we take i)0,Red = 16.8 ± 5.3 km s~^, 
from the average of the significant datapoints correspond- 
ing to the forbidden lines (Ha points are excluded because 
of the possible contamination with light coming from the 
star and accretion processes, and larger projection effects). 
From these values one obtains ro,Rod = 1-58 ± 0.36 AU. 

In the blue lobe the forbidden lines are quite faint in 
the 5*1 and 5*7 slits, but the needed velocity shifts can be 
extracted from the slit pair S6 — S2. Taking the average of 
the significant values measured within O'.'S from the star 
one obtains v^^biuc — 14.9 ± 4.5 km s~^ at 20 AU from 
the axis. With Wp.eiuo ^ 260 ± 10 km s^^, the origin of 
this portion of the flow appears to be located at r-Q.Biuc — 
0.44 ± 0.1 AU. The values for both the red and blue lobes 
are in good agreement with analogous determinations by 
Coffey et al. I)2004f) . We note that according to Pesenti et 
al. (2004) the toroidal velocity derived from measurements 
of the full velocity profile (as in our case) at robs = 20 AU 
from the axis may be underestimated by a maximum of 
15% because of projection effects. A 15% correction to the 
toroidal velocity for the blue lobe, however, is smaller than 
our accuracy and will not be considered in the following. 

It should also be remembered that the above calcula- 
tion strictly gives only the footpoint of the flow surface for 
which the rotational velocity could be measured, and not 
the outer radius of the whole ejection region. For exam- 
ple, Takami et al. H2004fl report about the discovery of a 
cold and slow wind component emitting in H2 lines, that 
surrounds the base of the optical jet from DG Tau. Such a 
component is probably anchored at larger footpoint radii 
than the optical component. Thus, the derived footpoint 
values may be conservatively considered as lower limits to 
the true extent of the launching region. 

4.2. Jet/disk angular momentum balance 

In the magneto-centrifugal approach, the angular momen- 
tum balance in all space can be written in conservative 
form as (Casse & Ferreira 120001 '): 

V • (^prv^vp - - rTv) = 0, (2) 

where T is the 'viscous' stress tensor inside the disk 
(Shakura & Sunvaev I1973|) . and use has been made of 
the relation V • rTv = r(V • Tv). Viscosity can transfer 
the excess angular momentum radially from the inner to 
the outer region of the disk, but it is very difficult to quan- 
tify its importance (see the discussion in Casse & Ferreira 
I2000|l . Our rotation measurements, however, allow us in 
principle to give rough estimates of the various contribu- 
tions as follows. 

We integrate Eq. El over the volume delimited laterally 
by the magnetic surfaces passing through ro^n and ro.cxt, 
and above and below the disk by two surfaces, Sucd and 
'5'biuc, parallel to the disk plane and located in the ob- 
served regions at ~ ± 0'.'4 AU from the source. We are 
thus led to calculate four surface integrals. By definition 
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Bp • n = along the magnetic surfaces in all space, and 
in the jet Tv = (ideal conditions), and Vp || Bp. Thus 
the integral over the lateral surfaces reduces to the fluxes 
through the annuli at ro^n and ro,ext- Since in the disk 
the poloidal velocity is equal to the accretion velocity, the 
sum of the kinetic component of the flux through the an- 
nuli under consideration is 

Min rosn WK,in " M^xt ?'0,cxt WK,cxt = — ^disk, (3) 

where wk = {GM^/ro)'^/'^ is the Keplerian velocity, and 
Afin, Mcxt are the mass accretion rates in the disk through 
the magnetic surfaces located at ro,in, rp^oxt respectively 
(cf. Fig. 6). Thus Ldisk represents the angular momentum 
that has to be extracted from the disk zone under con- 
sideration per unit time to allow for mass accretion. The 
component of the integral over the same surfaces due to 
viscosity will be indicated by Lt, while the angular mo- 
mentum fluxes in the jet lobes, ijct,R and ijct,B result 
from the calculation of the integrals over S'biuc and 5Rod- 
In summary, one reduces to: 

idisk = + ijct3cd + ijct, Blue- (4) 

Our rotation measures can be used to give estimates of 
the angular momentum loss rate in the disk, Ldisk, and of 
the angular momentum extracted by the jet per unit time, 
ijct. A comparison then between these quantities can give 
us a rough idea of the role played by the outflow, and, in 
turn, by viscosity, in the angular momentum balance. 

To calculate Ldisk we note that mass conservation leads 

to 

Moxt = Afjct,Rcd + Mjet,Bluc + A^in- (5) 

Referring to Fig.|Hl we assume that all the mass flowing 
through the annulus at rp^n is deposited onto the star. 
Thus a good estimate for the inner mass flux is Min = 
1.6 10~^ Mq yr~^, that was derived by Hartigan et al. 
( 1995|) from measurements of the veiling of photospheric 
lines caused by the accretion shock. The mass flux at the 
jet base in the two lobes, Afjet,Rod and Mjct,Biue, has been 
estimated in Paper I to be 0.6 and 1.1 10~^ M© yr~^, 
respectively. We thus obtain A/cxt = 1-77 10"^ M© yr'^ 
Taking conservatively ro^cxt = 1.58 ±0.36 AU, and ro^n > 
0.03 AU, one has 

idisk < 6.6 ± 1.5 10"^M©yr"i AUkms-^ (6) 

We will now attempt an estimate of ijct, Red and 
ijct, Blue- Since in the jet pvp = fc(a)Bp, this term can 
be written, for each lobe, as: 

ijct = {rv^ - pvp • n dS. (7) 

The MHD wind theory prescribes that the quantity in 
parentheses, which represents the total (kinetic plus mag- 
netic) specific angular momentum, is also constant along 
the magnetic surfaces : 

Zjct(a) = {rv^ - rB^/ink) = r\flo (8) 
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where = vx/ro is the angular speed of the disk at 
the base of the considered surface (see Konigl & Pudritz 
I2UUU|I . Equation [7| can thus be transformed to 



L 



'jet 



roUKpVp • n dS. 



(9) 



The quantity ^ represents the 'magnetic lever arm' of 
the flow, and it varies between 2 and 5 for typical flow 
parameters (see, e.g., Konigl & Pudritz f2U(JU[ Anderson et 
al. 120031 Pesenti et al. l2004|l . Assuming then that the flow 
surfaces are well-behaved and regularly nested, an upper 
limit to the angular momentum flux transported by the jet 
can be estimated from the measured jet mass flux rates, 
the derived values of vq and tjk at the outer footpoints 
(reported in Tables 13 and ^J, and assuming a value for the 
magnetic lever arm. For example with (r^i/ro)^ = 10 one 
obtains |Ljot,Rod| < 2.6 ± 0.8 10"^ Mq yr'^ AU km s"^ 
and |ijct,Biuc| < 2.5±1.2 10"^ M© yr~i AU km s^^. When 
summed, this amounts already to about 80% of the excess 
angular momentum, Ldisk, to be extracted from the por- 
tion of the disk from which the observed flow originates. 

An interesting lower limit to the angular momentum 
transported by the jet can instead be given for the self- 
similar disk wind model of Ferreira & Pelletier H1993|l . 
Before deriving this we note that the 3-D geometry of the 
B and v vectors in the magneto-centrifugal scenario is 
such that if the system rotates anticlockwise looking down 
the +z axis {v^,vk > 0), then B,j,Bp/pVp = B^/k{a) < 
in all space (i.e., in both the z > and z < domains), 
while the opposite holds for clockwise rotation. Therefore, 
for every possible configuration of the system the kinetic 
and magnetic terms in the specific angular momentum 
have opposite signs, or they add up in absolute value. 
Thus, for each lobe, a lower limit to |ijet| is set assuming 
that the magnetic contribution is negligible. The integral 
in Equation 13 can be computed exactly for the model of 
Ferreira & Pelletier l)1993|l . For this model v^,Vz (x r~^^^ 
and p oc r^^~^/'^^ , where ^ is the ejection index of the jet 
(see next section). In this case Eq. Oreads: 



1/2 
obs 



(C+1/2) _ (5+1/2) 
obs in 



e + i/2 



Mi 



jot, 



(10) 



' obs 



where the integral is calculated between the radius of the 
axial inner hole of the jet rin, that is about 5 AU at 0'.'4 
from the star (Cabrit et al. I1999|l and the radius of the 
observed external layer, robs, which is 30 AU and 20 AU 
for the red and blue lobes, respectively. With ^ = 0.040, a 
value similar to the one that will be derived for this jet in 
the next section, we obtain |Ljct,Rcd| > 2.4±0.8 10^^ Mq 
yr-i AU km s'^ and |ijct,Biuc| > 2.0 ±0.9 10"^ Mg yi'^ 
AU km s""'^. Again summing these quantities up the total 
amounts to about two thirds of the angular momentum 
-^disk, that has to be extracted from the examined portion 
of the disk per unit time. 



The rough estimates above indicate that the mecha- 
nism generating the jet may indeed be capable of braking 
the disk with high efficiency. In the rest of this discussion 
we will assume that in the system under study all of the 
excess angular momentum is carried away by the wind, or, 
equivalently, that in the disk the magnetic torque is much 
larger than the viscous torque. 

4.3. Magnetic lever arm and ejection index 



Following the disk-wind formulation by Ferreira H1997I 
2002), in a jet/disk structure like the one considered, the 
mass accretion flux through the disk. Mace, varies with 
distance from the star as Mace = Const x rg. The ex- 
ponent ^, that measures the local ejection efficiency, is 
called 'ejection index', and regulates many of the physical 
properties of the wind. Under the assumption of dominant 
magnetic torque (see above), ^ is simply related to r^/ro 
by the expression: (rA/j'o)^ = A 1 -I- 1/2^. 

According to the disk-wind models 0.005 < C < 0-5, 
where the upper limit applies to self-similar solutions with 
sub-Alfvenic heating (warm solutions, see Casse & Ferreira 
2000). It is interesting to find limiting values for this pa- 
rameter in our case. 

A lower limit for ^ is set by the conservation of mass 
in the disk/jet system. The external mass flux can be ex- 
pressed through ^ and Min as 



Mcxt = Min(7-0xxt/?'0,in)^ 



(11) 



Since ro^n cannot be lower than the corotation radius 
Tc, a lower limit for ^ is given by: 



ln(l + (Mjet,R,cd + Mjct,Bluc)/Min) 

ln(ro,cxt/?'c) 



(12) 



Taking for ro,oxt the value of ro derived for the red lobe 
(1.58 AU) one obtains ^min = 0.025 ± 0.005. If instead we 
take ro,cxt = 0.44 AU, that is the lateral footpoint value 
for the blue lobe, one obtains ^min = 0.037 ± 0.005. The 
value of ^min may further differ if one considers that the 
full extent of the accretion/ejection region may be larger, 
due to the possible presence of a surrounding cold flow 
not visible at optical wavelengths (such as a II2 wind, cf. 
Takami et al. [2004). This would imply not only a larger 
''o.ext, but also a larger Mjct, as one would have to consider 
also the contribution of such flow to the mass outflow rate. 
Both these quantities, however, are still undetermined for 
the case under study. 

An upper limit for ^ follows from the conservation 
of angular momentum in the jet region. For dominant 
magnetic torque the constant of motion can be expressed 
through A or ^ in terms of the Keplerian rotation velocity 
at the footpoint: 



IW(a)| 



B^ 
Ank 

1 



= A VKro 



^) (GM, 



(13) 
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Now, as described in the previous paragraph, the kinetic 
and magnetic terms in I always have opposite signs, so 
the condition jB^j > imposes an upper limit for ^. For 
the red lobe, using tq = 1.58 AU, one finds ^max,Rod = 
0.046 ± 0.011, while for the blue lobe, with tq = 0.44 AU, 
one obtains ^max.Biue = 0.041 ± 0.010. 

In summary, for the red lobe of the RW Aur jet 0.025 
±0.005 < ^ < 0.046±0.011, or 12±3 < A < 21±4. At the 
same time, in the blue lobe one finds 0.037 ±0.005 < ^ < 
0.041 ± 0.011, or 13 ± 3 < A < 15 ± 4. The derived values 
for ^ are thus in the lower part of the range admitted by 
the model. Nevertheless, "cold" wind ejection (enthalpy 
negligible with respect to the kinetic energy at the foot- 
point) requires ^ ^ 0.01 (Casse & Ferreira |2UU0|I . Thus 
our results are consistent with the idea that some sort of 
heating is provided at the base of the jet (i. e. the wind is 
"warm" ) , for example by a hot disk corona. The fact that 
an extra-heating seems to be required at this location has 
also been noticed by Garcia et al. l|2001a|l from the anal- 
ysis of terminal poloidal velocities and jet total densities, 
and by Pesenti et al. I|2004|l from the analysis of rotation 
signatures in the DG Tau jet. If the heating is produced by 
an active disk, however, a certain amount of viscous dissi- 
pation must be allowed. Alternatively, the heating at the 
jet footpoints might be due to ambipolar diffusion and/or 
X-ray irradiation. See Garcia et al. I|2001b|l and Shang et 
al. H2002|l for a discussion of these processes respectively. 

4.4. Properties of the magnetic field 

A quantity that can be derived directly from our mea- 
surements is the ratio B(f,/Bp between the toroidal and 
poloidal components of the magnetic field vector at the 
location of the observations. This quantity indicates how 
much the lines of force are wrapped on a given magnetic 
surface. In the magneto-centrifugal scenario, the coUima- 
tion of the jet is thought to arise from the hoop stress 
generated by the increase of the toroidal component of the 
magnetic field after the gas has passed through the Alfven 
surface (see, e.g. Konigl & Pudritz ,2000j . It is thus inter- 
esting to obtain from the observations an indication of the 
magnetic field configuration in the region just above the 
coUimation zone. 

The ratio Bff,/Bp can be derived using a further con- 
servation law of general disk-wind theory (see Konigl & 
Pudritz 2000, Anderson et al- HUU^ : 

[v, - ^) /r no, (14) 

where fio is the disk angular velocity at the footpoint. 
Introducing the values observed for v^, Vp at r = robs in 
the two lobes, and using the corresponding vq to calculate 
Qo, we find B^/Bp = 3.8 ± 1.1 in the red lobe, at 30 AU 
from the axis and about 80 AU from the disk. In the blue 
lobe, one finds B^/Bp = -8.9 ± 2.7 at 20 AU from the 
axis and at about 100 AU above the disk (the signs refer 
to a cylindrical coordinate system with the positive z axis 
oriented along the blue lobe). 



11 

If the poloidal magnetic field is almost entirely pro- 
jected along the z direction, these ratios would correspond 
to an inclination angle of the magnetic line of force with 
respect to the disk plane of about 15 degrees for the red 
lobe, and only 6-7 degrees for the blue lobe. The dom- 
inance of the toroidal component of the field appears to 
confirm that the flow is magnetically coUimated, as pre- 
scribed by the magneto-centrifugal models. Also, the well- 
known asymmetry between the lobes of this jet is reflected 
in the magnetic configuration, as the field turns out to be 
more tightly wrapped in the blue lobe. 



5. Conclusions 

In this paper we have described a new combined observa- 
tional/theoretical study of the rotation properties of the 
bipolar jet from RW Aur. 

We have found rotation signatures in a set of spectra 
taken with HST/STIS with multiple slits oriented paral- 
lel to the flow axis, using techniques similar to the ones 
adopted in our first rotation study (of the DG Tauri jet, 
Bacciotti et al. lMT^ . We analyse the first 300 AU of the 
jet (1'.'5 at the distance of RW Aur), applying an updated 
version of the correction routines for uneven slit illumina- 
tion. We find that both lobes rotate in the same direction 
(i.e. with different helicities), with toroidal velocities in 
the range 5-30 km s~^ at 20 and 30 AU from the sym- 
metry axis in the blue and red lobes, respectively. The 
sense of rotation is anti-clockwise looking from the tip of 
the blue lobe (P.A. 130° North to East) down to the star. 
These results are confirmed by other HST/STIS observa- 
tions from our group (with the slit placed transverse to 
the jet axis) presented in Coffey et al. 1)2004 1) . Rotation 
is more pronounced in the data from the outermost slit 
pair, as expected because of projection effects along the 
line of sight (Pesenti et al. "20041. Also, rotation signatures 
are more evident in [01] and [Nil] lines than in Ha and 
[SII] lines. We interpret this result as due to the fact that 
[01] and [Nil] emission traces regions closer to the jet axis 
than Ha and [SII] , and hence at a given position of the slit 
the jet is observed closer to its border in [01] and [Nil], 
which reduces the confusion due to projection. 

The observed rotation favours widely known magneto- 
centrifugal models for the jet generation. Following the 
formulation in Anderson et al. f2003i , the derived toroidal 
and poloidal velocities imply that the fiow surfaces of the 
redshifted lobe observed at 30 AU from the axis are rooted 
in the disk at about 1.6 AU from the star. The blueshifted 
flow observed at 20 AU from the axis arises instead from 
a region in the disk at about 0.4 - 0.5 AU from the star. 

Making use of general principles underlying the disk- 
wind models, we have also derived other parameters use- 
ful to constrain the properties of the RW Aur accre- 
tion/ejection structure. We have estimated upper limits 
for Ljot, the angular momentum transported by the vis- 
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Table 3. Measured physical quantities at the jet base 



Jet lobe 


robs 


a,b 

V , 

4> 












-'"jet 




(AU) 


(km s 




(km s" 




(Mo) 


(10~6 Ms yr-i) 


(10-« Mq yr-1) 


Redshifted 


30 ± 3 


16.8 ± 


5.3 


-130 ± 


10 






0.06 ± 0.01 














1.3 ± 0.1 


1.6 ± 0.1 




Blues hifted 


20 ± 3 


14.9 ± 


4.5 


260 ± 


10 






0.11 ± 0.01 



Average of the significant values obtained from the FELs between 0'.'2 and 0'.'6 (O'.'S) for the red (blue) lobe. 
^ The vector components are in a cylindrical coordinate system centred on the star and with positive z axis oriented along 

the blue lobe (South-East in Fig. 1). 
" From Woitas et al.. .2001. 

From Hartigan et al. 

From Paper I. 



Table 4. Derived physical quantities 



Jet lobe 


(AU) 


(km s-^) 


1 -^'jct ,min 1 

(Meyr-^AUkm/s) 


-'-^disk.max 

(Moyr-^AUkm/s) 


r 




B4,/Bp 


Red 
Blue 


1.58 ± 0.36 
0.44 ± 0.10 


27.0 ± 3.2 
51.2 ± 6.2 


2.4 ± 0.8 lO"'^ 
2.0 ± 0.9 lO"'^ 


6.6 ± 1.5 10"^ 


0.025 <i< 0.046 
0.037 <i< 0.041 


12 < A < 21 

13 < A < 15 


3.8 ± 1.1 
-8.9 ± 2.7 



^ Derived with the method in Anderson et al. 120031 
^ Keplerian velocity at the footpoints ro. 

Imposing that \B^\ = at the location of the observations and adopting the scalings in Ferreira & Pelletier 1993 (see text). 

Considering only the portion of the disk from where the outflow seen at optical wavelengths originates. 

Assuming dominant magnetic torque; accuracy 20 - 22 %. 
' At about 80 AU above the disk and 30 (20) AU from the axis for the red (blue) lobe. Coordinate system as in Tabled 



ible jet lobes, and a lower limit for the same quantity 
in the special case of the self-similar disk-wind model of 
Ferreira & Pelletier ((1993.) . We compare these values with 
the angular momentum, Ldisk, that the region of the disk 
from which the visible outflow originates has to lose per 
unit time in order to accrete at the observed rate. We 
conclude that the jet is capable of extracting a consistent 
fraction (two thirds or more) of the excess angular mo- 
mentum present in the disk. 

Assuming moreover that all the excess angular momen- 
tum is carried away by the jet, we have also estimated the 
magnetic lever arm (expressed by the ratio r^/rg between 
the Alfven and footpoint radii) for the self-similar disk 
wind model of Ferreira & Pelletier H199l-{|l . We found this 
quantity to be in the range 3.6 - 3.9 for the blueshited 
lobe, and in the range 3.5 - 4.6 in the redshifted lobe (ac- 
curacy 20 - 25 %). Alternatively, the value of the ejection 
index ^ — dln(Macc)/rf?' varies from 0.025 to 0.046 in the 
red lobe, and from 0.037 to 0.041 for the blue lobe (with 
the same accuracy). We caution that our determination 
of the magnetic lever arm (or ejection index) may be af- 



fected by the poor knowledge of the global extension of 
the launching region. In fact a component of the flow too 
cold to be visible at optical wavelengths may surround 
the observed jet. Nevertheless, the values determined for 
the optical component are in the range predicted by MHD 
models, and they also suggest that some heating is pro- 
vided externally at the base of the flow. The nature of 
such heating, however, remains to be identified. 

Finally, we have used our rotation measurement to 
derive information about the spatial configuration of the 
magnetic field in the examined region. In particular, us- 
ing well-known conservation laws of the MHD disk wind 
theory we have derived the ratio B^/Bp between the 
toroidal and poloidal components of the magnetic field 
at the observed locations in both lobes of the bipolar jet. 
We obtained B^/Bp = 3.8 ± 1.1 for the red lobe at 30 
AU from the axis and about 80 AU from the disk and 
B^/Bp = -8.9 ± 2.7 for the blue lobe, at 20 AU from 
the axis and 100 AU from the disk (in cylindrical coordi- 
nates, with positive z along the blue lobe). The toroidal 
component of the magnetic field appears thus to be dom- 
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inant, as expected for a magnetically coUimated flow. In 
addition, the field seems to be more tightly wrapped in 
the blue lobe, reflecting the well-known (but unexplained) 
asymmetries between the two lobes of this jet. 

In summary, our observations and subsequent analysis 
appear to conflrm once more the magneto-centrifugal sce- 
nario for the launching of YSO jets. To prove conclusively 
that jets solve the "angular momentum" problem in star 
formation, however, will require further detailed studies 
of a larger number of jets, and possibly at higher spectral 
and spatial resolution. 
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